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Crack length 
Thickness 
Young's modulus 

A measure of strain hardening in the assumed bilinear stress-strain 
behavior of the material vide eqiiation (1) 

Compliance function^ £vB/P 
K calibration factor » K/PB^^ 

Stress intensity factor 

Value of stress intensity at which the crack starts growth 
ASTM £399 valid plane strain fracture toughness 

is the maximum value of stress intensity used during the final stages 
of fatigue cracking 

Plasticity corrected stress intensity factor value corresponding to 
a given crack extension as in the R-curve determination 
Stress intensity factor corresponding to 2% crack extension as mea- 
sured by the secant technique 

Stress Intensity factor corrected for plasticity as per equation (4) 
Yielding constraint in K-R relationship in equation (2b) 

Yielding constraint in K-Rg£f relationship in equation (2a) 

Load 

Distance from the crack tip along x direction 
Plastic zone size 
Nominal plastic zone size 
Width 

Critical width at which the crack initiates 
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Normal stress in y direction 

Uniaxial yield strength determined according to .02 offset procedure 
Normal strain in y direction 
Yield strain and is equal to Oy/E 
Total displacement 

Elastic component of the displacement 

Total deviation fr<» the linear displacement 

Deviation from linear displacement due to the growth of plas-* 

tic zone 

Deviation from linear displacement due to the growth of crack 
Total crack extension 
Physical crack extension 

Crack extension equivalent to the growth of plastic zone 
and is equal to Rg^^/2 



A NEW BASIS FOR THE DETERMINATION OF FRACTURE TOUGHNESS 


S. Banerjee* 

Ames Research Center 

INTRODUCTION 

Plane strain fracture toughness » ^ property of a material and 

is expected to be independent of the size, configuration, and loading of the 
specimen or the structure. A procedure for the determination of is 

described in ASIM E399 In this procedure, the start of crack growth 

is identified by the deviation from linearity in the load-displacement test 
record. This approach is simple and reasonable. But plasticity preceding 
or accompanying the crack growth is not considered in a satisfac- 
tory manner which can cause Kq or, in some cases, to depend on the 

specimen width. 

A new approach is proposed here which gives a width-independent 
value* In the new approach, is defined as the value of the stress 

intensity factor, K, at which the crack starts physical extension. The 
approach is based on the results of approximate analyses and supporting 
experimental data* At first, the existing approaches to the determination 
of are examined. Then, an analysis is presented which demonstrates 

that the growth of the plastic zone and the constraint in a compact tension 
specimen can depend significantly on the specimen width. Afterward, a 
typical R-curve [ 2 ] is analyzed to evaluate the simultaneous contribution 

of crack growth and crack tip plasticity to the observed deviation from 

^Senior NRC Resident Research Association on leave from the India 

Institute of Technology, Bombay. 

^The italic numbers in brackets refer to the list of references 


appended to this paper. 



linear elastic behavior. A coodi>inatloa of these two analyses produces a 
very simple procedure for the determination of fracture toughness. The 
final discussion includes a comparison of the analytically predicted and 
experimental results* a description of the proposed approach, and the 
implications of the results obtained. 

EXISTING APPROACHES TO THE DETERMINAT1(»I OF 

xo 

In ASTM E399 [I], tte load, P^, corresponding to the start of crack 
growth is identified by the 5Z secant technique. The 5Z secant is supposed 
to correspond to an equivalent crack extension of 2%. 

According to ASTM E399, three possible types of load-displacement test 
records are obtained during the testing of the precracked specimens. A 
material with low toughness and high yield strength in which the crack 
extension occurs abruptly gives what is termed a Type-III load-displacement 
test record [I]. The value of calculated from Pq in such a case is 

largely Independent of the size of the specimen. The Type-II load- 
displacement test record is similar to the Type III in this respect. How- 
ever, in materials with somewhat higher toughness, the crack and plastic 
zone grow simultaneously. This type of behavior gives a Type-I load- 
displacement test record. According to E399, in such a case, the value of 
Kp computed from Pp is considered to be a valid measurement if the 

following two conditions are satisfied: (a) the specimen thickness, 

B 2 2.5 (K»/a^)^ where a is the yield strength; and (b) P /P_ < l.I 
X ^ I max y 

where P is the maxlmim load encountered during the loading of a pre- 
max 

cracked body. For clarity, the obtained in such a manner will be 



referred to as [K^„l . Recent experiments [3-5] have shown that K- and, 

IC 399 Q 

in some instances, the valid [K_^] depend on specimen width W. 

i-C 3 99 

ASTM E399 implicitly assumes that the plastic zone size depends only 
on specimen thickness and is Independent of width. On the other hand, the 
width dependence of K_ or [K _] values is observed in the compact ten- 

0 It> 3 99 

slon specimens which satisfies the ASTM E399 requlranent that the thickness 

>2.5 (Kq/o^)^. Further, the experimental data indicate that the thickness 

has only a small effect on values obtained from the compact tension 

and three-point bend specimens. It thus appears that a large thickness 

cannot guarantee a limited plasticity during the start of crack growth unless 

the width is correspondingly increased. 

The second condition above, concerning P /P., is specified to 

max 

guarantee that the deviation from linearity is produced primarily by 

crack extension. The condition that P /?_ < 1.1, can be achieved 

max V 

only when the R-curve [ 2 ] is flat; where the R-curve is a plot of stress 
intensity factor corrected for plasticity vetsus the corresponding physi- 
cal crack extension. For a given material, a load-displacement test 

record which satisfies the requirement that P /P^ S, 1.1, can be obtained 

max Q 

in a very limited combination of specimen size and configuration. A speci- 
men is tested in a laboratory to simulate the behavior of a structure in 
service. Yet the form of load-displacement test record of the specimen 
where Pmax/^Q ^ correspond to the form observed in a real 

engineering structure. Therefore, [Kt/.] of a material measured in a 

it 39 9 

specimen where ^niax^^Q ^ meaningful in predicting the 

fracture behavior of a structure where Pjnax^^Q ^ 
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The problem of the width dependence of Kq values has been 
examined by others. Kaufman [5] has proposed a relaxation of the 
/P. requiresffint and the use of specimens with thickness/width 

IBaX ^ 

ratios of less than 0.5. Newman [101 has suggested a two-parameter 
approach to predict toughness In specimens with varying widths. This 
approach is empirical and is concerned with the maximum load at failure 
and therefore relates to the final propagation rather than the start o;-! 
the crack growth. Munz [J] has proposed an approach idiere the toughness 
is based on the actual point of crack extension. He has recommended a 
variable secant technique where the secant value is adjusted according 
to specimen width. 

PLASTICITY, CONSTRAINT, AND WIDTH OF A COMPACT TENSION SPECIMEN 

Recently, it has been shown through a simple and approximate analy- 
sis [ii-I3] that the width and crack length of single-edge notch speci- 
mens at a given value of K significantly influence the nominal plastic 
zone size, the yielding constraint, and the crack opening displacement 
(COD). These results are in qualitative agreement with the recently 
reported experimental observation that the value of Kq increases with 
increasing width [3-d]. The width dependence of Kq suggests that the 
growth of the plastic zone ahead of the crack tip may be dependent on the 
specimen width. On the other hand, the two-dimensional finite element 
analyses do not give entirely satisfactory results. For Instance, the 
plastic zone size can be considered equivalent to a virtual crack exten- 
sion. Accordingly, the displacement of the specimen depends on the 
plastic zone size. But the finite lement analyses give width-dependent 



plasticity-induced displacement values that are substantially less than 
the experimentally determined ones [3] particularly when the plastic zones 
are small. Thus, the reported two-dimensional analyses do not give a sat- 
isfactory representation of the growth of small plastic zones. In fact, 
it must be recognized that the crack tip deformation in the fracture tough- 
ness specimen is a three-dimensional elastic-plastic boundary value problem. 
Since such a problem is not satisfactorily resolved as yet. the simple and 
approximate analysis reported earlier [11-13] is used to evaluate the 
effect of specimen width on the displacement, constraint and the Kq value. 
What follows in this section is a brief outline of this analysis and a com- 
parison of the experimental and the analytical results. The good agree- 
ment between the experimental and the analytical results validate the analy- 
sis reported here. 

The loading of a compact tension specimen can be satisfactorily rep- 
resented by the combination of an axial force applied at the midpoint of 
the ligament (W-a) and a bending moment [14], An (r)~^^^-type strain dis- 
tribution at distances close to the crack tip is assumed, where r is 
the distance from the crack tip. The assumption follows from experimental 
observations [15-17]. The linear (-r) type distribution arises because 
of the bending of the compact tension specimen [id]. Experimental obser- 
vations confirm that such a composite strain distribution exists ahead of 
the crack [19], The composite strain distribution is schematically pre- 
sented in Fig. 1. a plot of local strain Cyy against the distance X. 
measL -sd from the tip of a crack of length a. The strain. Cyy. is related 
to the stress. Oyy» through an assumed bilinear stress-strain behavior of 
the material which can be represented by 
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a ■ e E for e < e_ 

yy yy yy y 


and 


( 1 ) 


e E + (e - e_)E’ • + (e - Gv)®’ fo*' 

yy y yy x y ' yy 


yy 


where E* is a measure of the strain hardening of the material. In the 
subsequent restilts calculated and reported here, E* is assumed to be equal 
to E/150. It may, however, be noted that the results correspond to a 
limited amount of crack tip yielding, and therefore are not significantly 
Influenced even if E* is assuaed to have a value 3 times less, that is 


E/450. 


As mentioned above, the analytical results reported here correspond 
to a limited amount of crack tip plasticity. Therefore, the plastic zone 
and the strains at the crack tip are rather small. Thus, the conversion 
of the strain into stress through equation (1) is a reasonable procedure. 

The bending moment and the axial force can then be written in terms of 
the 0 stress. 

yy 

Based on the composite strain distribution (Fig. 1), load and moment 
balance equations were formulated. These equations were solved using 
Brown's method [20] • From this, the strain reversal point, Xa, and the 
plastic zone size, R > (Xz - Xi), were obtained. The loadline and crack 
mouth opening displacements of a compact tension specimen cooq>uted accord- 
ing to the above procedure agree well with the experimentally observed 
values [12], 

* 

The load-moment balance equations were solved for progressively 
increasing loads for a given specimen size. The calculations were performed 
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for different widths at a/w « 0.4> O.S, and 0.6 to evaluate the effect 
of width and crack length on the value of R at different loads. 

As shown in Fig. 1, the nominal pl.^stic zone size, R, is defined by 
the point ahead of the crack tip at which “ ^Y* 

yield strain and is equal to a^/E. Therefore, at r » R, * o^. 

Figure 2 is a plot of the strain, versus r, the distance from the 

crack tip in the plane of the crack. V/hen the load is small, the plastic 
zone size has the shape of a pair of horns as schematically presented in 
Fig. 2. The analytical results show that in such a situation the actual 
plastic zone size is approximately equal to R, the distance from the 
crack tip at which 121 ]. Correspondingly, that is in the plane 

strain condition, the displacement and the change in compliance is pro- 
duced by an effective plastic zone size, Eg££t which is less than R [22] 
as shown in Fig. 2. Both R and Rg££ depend on the state of stress or 
the yielding constraint at the crack tip. 

According to McClintock and Irwin [ 22 ], R^££ is related to K 
through a relationship 

Reff “ (K/ay)^/Mir (2a) 

where M is a measure of the yielding constraint. Likewise, R is 
related to K through a similar relationship 

R " (2b) 

where m is also a measure of the yielding constraint. It is assumed 
in this paper that 

R ff - R/m (3) 
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Substitution of Eqs. (2a) and (2b) in (3) gives M • Thus, m and M 

can be calculated if K and R are known. 

The strain, e , at a distance, r > r ahead of the crack tip can 
yy eff 

be obtained from Eq. (3) and the 1//F dependence of the strain as 

shown in Figs. 1 and 2. The 1//T dependence of the strain gives 


(e ) 


(e ) 

yy 


^eff / R Y 

..a ■ 


Accordingly, 


(e ) ” /m(e ) » /me„ 

!'>’ r-R „ r.R ^ 

eff 


This result is shown schematically in Fig. 2. 

The assumption of Eq. (3) is reasonable and is consistent with the 
following observation. At a stress Intensity value corresponding to 
equivalent to 2% crack extension, the value of M calculated according to 
the analysis reported above is around 3. At higher stress Intensities, 
the plastic zone size Is large and the shape of the plastic zone tends to 
be circular. The state of stress In this case corresponds to the so-called 
plane stress case. Correspondingly, * R and M « m ■ 1. These 

observations are In general agreement with the proposed values of M for 
plain strain and plane stress conditions [22]. In addition, at r ■ R^££, 
Cyy ■ /m Cy and correspondingly / m o^. This result Is consistent 

with the analytical results on elastic-plastic stress distribution and the 
nominal plastic zone size [21], 

The results from our analysis are next examined through two different 

relationships: K*la„/v as a function of AV /V and K*/a /)7 as a function 

Y p e y 
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of c.. The former is concerned with the growth of the plastic zone and 
the latter with the change in the yielding constraint. K* is the stress 
intensity factor, K, corrected for plasticity according to equation (4) 
discussed presently. 


K*/o^/w and the Growth of the Plastic Zone 

In the absence of crack growth, the growth of the plastic zone 
directly relates to the deviation from the linear elastic behavior. 
According to Irwin's concept of virtual crack extension [23], the effec- 
tive crack length Is given by 

®eff *^eff^^ 

where a is the physical crack length. The values of <^sn be 

obtained from the solution of the load and the moment balance equation. 

The effective crack length is then obtained from equation (4) and is 

used to calculate K* which is called the plasticity corrected K value. 

The stress intensity factor, K., is calculated from the physical crack 

length, a. It may be noted that the difference between K and K* is 

rather small since the ir ^.Imum value of R that could be obtained from 

eff 

the analysis is small. 

Values of R ,, can be used to calculate the deviation from 
eff 

linearity through the standard compliance relationship. The deviation 

from linearity is expressed in terms of it*® terms 

and dV are the linear elastic displacement and the deviation from 
P 

linearity produced by the growth of the plastic zone. The terms V, 

4V, and AV^ are the total displacement, the total deviation from lin- 
earity and the deviation from linearity produced by the growth of the 
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crack, respectively. All these terms are schematically represented 
in Fig. 3. 

An examination of the results obtained from the analysis outlined 
in Figs. 1 and 2 shows that the parameter K*/o^/kT depends on the value 
of AV /V as shown in Fig. 4. 

The AV^/V value, for a plastic zone size equivalent to 2Z crack 
extension, depends on a/W. This can be independently co^mted from 
the standard compliance relationship [24] z 

(BB/P)V « f (5> 

%ihere E is Young*s modulus, B is thickness, P is load, and i Is 
a function of a/W only. Differentiation of Eq. (5) gives 


^ AV 
P “ p 


df * 

— Aa 
da p 


d(a/w> W 


df 


Aa_ 


( 6 ) 


where Aa is the virtual crack extension due to the growth of the plas- 
P 


tic zone. As outlined in Eq. (4) and Fig. 3, Aa 


*eff/2* ®^^^*** 


Eq. (6) by (5) and obtain 


AV 

V 


lllfE 

f w 


r« Aa _ 
f ’ pa 


(7> 


e 

Table 1 lists AV /V values for Aa /a - 0.02 as calculated from 
p e p 

Eq. (7) or the different a/W values. The corresponding values of 
K*/o^^ are obtained from Fig. 4 and are also reported in Table 1. As 
shown in the tabl2, for equivalent to 2Z crack extension, the cor- 

responding value of K*/ay»^ is approximately equal to 0.5 for the a/W 
ratios investigated. The K* value obtained from K*Oyi^ ■ 0.5, corresponds 
to the value obtained by the secant technique of ASTM E399 and therefore 
K* = Kq. Thus, the analysis shows that • 0.25. This result is 

later compared with the experimental values obtained for different widths. 
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K*/o^i^ and the Yleldhig Constraint 

The values of K* and R obtained fron the results of the analysis 
outlined in Figs. 1 and 2 can be used to detemtne values of • and M 
froB Eqs. (1) and (2). As we have seen, both ■ and N are relative 
■easures of the yielding constraint. 

Variations in the paraseter, with b for various values 

of a/U are shown in Fig. 5, which shows that as II is Increased, b 
will increase. An increase in the value of b increases the crack tip 
local stresses and this, in turn, will cause a typical plane strain-type 
fracture. FroB Fig. 5 for a/W » 0.5, b » 1.9 at K/o^#^l * 0.5. The 
corresponding value of M will be about 3.6. 


Calculated and E]q>eriBental Values of the Paraseter 




It was shown that » 0.25 corresponds to a plastic zone 

size equivalrat to 22 crack extension and signifies a value of the yield- 
ing constraint which is constant at a given a/W. Figure 6 coBpares the 
experlnental and the calculated relationship, (Rq/Oy)^, versus W for a 
variety of Baterials: low-strength steels [J]^ and aluBlmai alloys 

[J, f]. nie Kq values reported are obtained under conditions where crack 
growth is absent. The scatter of the data in the figure aay be caused by 
the usual errors in the deteraination of Kq. 

Errors In the experlBentally deterclned Kq value can arise due to 
the following reasons: (a) the friction at the loading pins and at the 

supports for the clip gauge affects tlie load-displacement test record; 


^S. Banerjee, unpublished work. 
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(b) the secant value used to determine depends on a/W and should 

be calculated from Eq. (7). Instead, quite often, a 5Z secant as recom- 
mended in ASTM E399, is used for all a/W ratios; (c) Inadequate accuracy 
of the load-displacement test record introduces error in the secant mea- 
surements; and (d) the maximum stress intensity value used for fatigue 
cracking, K^(max) [ 7 ], influences the Kq value determined. 

The effect of K^(max), the maximum stress intensity value used in the 
final stages of fatigue cracking on R^, was investigated and is reported 
in Table 2. The values were measured according to the general pro- 

cedure outlined in ASTM E399 [ 7 ]. The friction at the loading pins was 
minimized by the use of flat bottom clevis and that at the clip gauge 
support was reduced by the use of razor blades. The load and secant 
value used to obtain were determined from Eq. (7), and the load and 

the displacenents were recorded with an accuracy of ±0.25Z. The 
value %ias determined with an estimated accuracy of ±2Z* 

As shown in Table 2, R^ values can be significantly increased due to 
the crack closure that occurs when the plastic zone formed at the crack 
tip during fatigue cracking is unloaded. Obviously, the R^(max) value 
should be low, preferably below < 0.22 to ensure the 

absence of the effect of R^(max) on the determined Rq values. 

Considering the potential sources of error in the determination of 


Rq, the agreement between the experimental and calculated values in 
Fig. 6 is good. The agreement implies that Rq/o^*^ Is a constant and 
is equal to 0.5. For convenience in future discussions, the calculated 


line in Fig. 6 is referred to as "line A." Thus, in the absence of crack 


growth. Rq 


values measured correctly should fall on line A. 


values 
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■easured in the presence of crack growth will be located to the right of 
line A. 

It is obvious froB Figs. 4, S» and 6 that, as W increases at a 
given a/V, Kq will increase, but M will remain constant. Finally, 
at a critical W • W^, Kg will approach the stress intensity value at 
which the crack starts extension, This is termed as the fracture 

toughness of the material. Std>stitution of V * and 
relationship Kq^/Oy^W - 0.25 gives • 4Kjg/o^*. At « > Wg. Kq > K^^,, 
and the crack growth starts. At W < W^, only the plastic zone grows as 
W increases. This can be represented by line A in Fig. 6. 

R-C08VE APPROACH AT H > 

The contribution of a limited amount of crack growth to the devis'- 
tion from linearity is evaluated from the analysis of a typical R-curve 
data [2], The contribution is evaluated in two different cases: 

(a) crack growth in the absence of plasticity (where the ct. < extends 
abruptly such as in a Type III or II load-displacement test recoru [I]>, 
and (b) crack growth together with the growth of the plastic zone (the 
case encountered in Type-I load displacement test record [i]). 

Crack Growth In the Absence of Plasticity 

Figure 7 is a schematic representation of a typical R-curve for a 
specimen with width W > U^, where K|^ is the stress intensity value 
corrected for plasticity for a given value of crack extension, Aa. The 
R-curve is a measure of the resistance of a material to crack growth. The 
value for zero crack extension is equal to the fracture toughness 


Kq-Kjc 
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defined in cMs paper. It is assumed that the R-curve for a limited 
amount of crack extension, that is, da/a ^ 0.02, can be represented by 



oc 


( 8 ) 


Tile assiffliption in Eq. (8) has the following justifications: (a) the 

It. rm of the equation is dimensionally homogeneous; (b) the R-curve can 
al"-o be represented as a plot of the parameter C versus Aa, where 
0 For small amounts of crack extension, the relationship in Eq. (8) 

is a good approximation of the typical R-curve data represented in terms 
of G [25]; and (c) the term on the left-hand side of Eq. (8) is a measure 
of the increments in the energy dissipated after the crack growth starts 
and therefore should be related to the amount of crack extension, Aa. 

* Kf.. Thus, Eq. (8) can be written as 


At 1 % crack extension. 



a Aa 


( 9 : 


Dividing both sides by W^, shown to be a property of the material, we 
obtain 


K 2 _ v2 

0 Aa 


<7 

«Y "c 


W. 


( 10 ; 


At Aa/a > 0 02 and a/W ■ 0.5, Aa > 0.01 W. In addition, as explained 
earlier, tne crack extension will start at a point on line A (Fig. 6). 
Thernrore, * 0.25 substitution of these two condi- 

■»ns in Eq. (10) gives the following proportionality: 


U 



How* 


w 

V "c 

Differentiation of the above equation gives d(KQ^/K|^)/d(HAf^)« 
ever, at Aa/a " 0.02, > K, and Aa « 0.01 H since a/tf • 0.5. 

Therefore, 


a(w/w^> 


0.01 



Aa/»^).02 


(ii; 


The ten on the right-hand side of Bq. (11) can be evaluated, at aiqr given 
values of and and a/U ratio, fro* F, the standard K calibration 
factor for a given speciaen geoaetry and a/V value, that Is, 



Aa/a*0.02 


1 

J 


»<^>1 
,2 >(4a/W)| 


Aa/a*K).02 


(12) 


where F >■ KBi^P and Is a function of a/H only [2d] and 

Fda - F((a + Aa)AlJ 


The K calibration factor, F, as a function of a/H Is known [2d]. 
Therefore, tl» ten on the right-hand side In Eq. (12) can be evalu- 
ated nuaerically. The nvaerlcal calculations show that the ten, 

1/F^. 9(F^)/9(Aa/H) does not vary much, over a considerable range of 


crack growth, that Is, at 0 < Aa/a < 0.05. At a/H “ 0.5 and 
Aa/a " 0.02, this ten has a value 

1 »(^L) , ■ 

F2 3(Aa/H) 

Substituting Eqs. (13) and (12) in (11) gives 

3(K^2/iq ) 


(13) 


(lA) 
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The 2Z crack extension in the absence of crack tip plasticity can there- 
fore be represented by a straight line with a slope of 0.064 in the 
(Kp/Kjp)^ versus W/W^ space. This line is termed "line C" in Fig. 8. 

Line A from Fig. 6 is also shown in this figure. The slope of line A is 
1 in Fig. 8 since » (4KQ^)/(a^^W^). The intersection of 

lines A and C gives of the material in a situation where crack 

extension occurs abruptly without any plasticity. 

Two observations can be made from Fig. 8. First, the start of crack 
extension preceded by the growth of a plastic zone equivalent to 
Aa/a = 0.02 results in approximately a 16-tlmes decrease in slope in 
the (Kp/Kj^)^ versus W/W^ plot. Therefore, the point of crack initia- 
tion can be readily Identified if the data were represented in such or 
an equivalent plot. However, if a Type-III load-displacement test record 
[1] is obtained during the test, plasticity preceding crack growth is 
less than 4a^/a = 0.02; in such a case, the decrease in slope in the 
load-displacement test record as a result of crack initiation would be 
quite drastic and the point of crack initiation can be readily identified. 
Second, line C is almost parallel to the x axis which represents W. 

This indicates that, in the case where the crack extension occurs abruptly, 
an increase in W has little effect on K determined by the intersec- 
tion of the line C and any other line A with an arbitrary slope value 
less than 1, in Fig* 8* Line A with a slope value less than 1 indicates 
that plastic zone preceding crack growth is less than an amount equivalent 

to Aa /a « 0.02a And this can occur in specimens with W > W . The 
P c 

above discussion agrees with the observation that is either independent 
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of the width or only mildly dependent on it. in the case where a Type-lII 
load displacement test record is obtained. 

However, it must be noted that, during fracture toughness testing, 
quite often the crack grows together with the growth of plastic zone. 

We have to consider, therefore, the sioniltaneoue growth of the crack 
and the plastic zone. This is done presently. 


Crack Growth in the Presence of Plasticity 

In the case where the crack grows together with the plastic zone, 
the total deviation from linearity, AV/V, consists of two terms: 

(a) the term results from plastic zone growth; and 

(b) the term that results from crack growth. These terms were 

earlier defined in Fig. 3: 



(15) 


The relationship between AV^/V^ versus was shown earlier in 

Fig. 4. At Aa/a * 0.02, K « Kg; therefore the relationship can be 
approximately represented by 


AV K 2 

-7~ *0.2 — y - for a/w 


0.5 


(x6) 


The contribution of crack extension to the deviation from linearity is 
calculated from the standard compliance relationship in a manner 
described earlier in Eqs. (5), (6), and (7). Accordingly, 


AV 


11 ^ 

f W 
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At a/W "0,5, f/f can be evaluated from the reported compliance 
relationship EVB/P - f as a function of a/W [ 24 ]. At a/W ■ 0.5, 
f'/f ■ 5. Therefore, the above equation reduces to 


AV Aa 

— 5-jc 

Vg w 


(i:*) 


It was earlier indicated that the slope, as given in Eq. (13), does not 
change much with Aa/W. Besides, at Aa/a • 0.02, K ~ Kq; one can then 
write 




Aa^ K^2 ya(K2/o^2„) 

"vT “ 


K 2 


Oy2w/ a (Aa/W) 6.40 y2W 


at a/W *0.5 


(If) 


Substitute Eqs. (18), (17), and (16) into (15) and solve for 
One then gets 


Kq2/o^2 „ 0.051 W 


(19) 


Divide Eq. (19) by 0.25 W^ and obtain 



( 20 ) 


Thus, the combined effect of the crack growth and the plastic zone can 
be represented as straight lines with slopes equal to 0.051 in a 
Kq2/0y^ versus W plot and 0.204 in a Kg^/K^j, versus W/W^ plot. 
For convenience, these lines are referred to as line B. 


DISCUSSION 

Comparison of Experimental and Calculated Results 

The Kg or values have been determined in the CT specimens 
for different combinations of progressively Increasing widths and thick- 
nesses prepared from titanium alloys [ 4 ] and an aluminum alloy [d]. T'ne 
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experimental and the calculated values of the parameter are 

compared in Fig. 9. 

The experimental Kq or values fall around the conq>uted lines 

A and B. This has some Interesting implications. Some of the experi- 
mental data points lie to the left of line A. The reason for such 
scatter In the experimental data has been explained earlier during the 
discussion of Fig. 6. Generally, the agreement between the experimental 
data and the calculated line B is quite good. The reasonably limited 
scatter around the line A and the good agreement of the experimental data 
with line B is particularly striking when one considers that the experi- 
mental data were obtained with specimens whose thickness varies by almost 
one order with the B/(Kq/Oy)* value ranging approximately from 0.5 to 5. 
If the constraint were to depend on the thickness as assumed in ASTM E399, 
Kq should change significantly with thickness. 

It should be noted that the intersection of lines A and B gives a 
Kq value equal to of the material. 

The results of the analyses presented in this paper predicts that, 
in a (Kq/Kj^) 2 versus W/W^ space, the experimental data point for 
all materials should fall around a single B line. Accordingly, the 
experimental data points located close to the B lines in Fig. 9 are 
replotted in Fig. 10. The experimental data, obtained from materials 
with different K^^ and a^, give results that fall around a single 
line B. The agreement indicates that the analyses of the width depen- 
dence of toughness and the assumptions made in the present paper are 
essentially correct. As pointed out earlier, this agreement is all the 
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more notable when one considers that the data obtained from the specimens 
have thicknesses that vary widely. 

A New Approach to the Determination of 

As Indicated earlier, in the present approach is defined as 

the oLress intensity value at which the crack extension starts. Since 
the start of crack extension being considered in this paper is that which 
satisfies ASTM E399 requirements, this will be preceded by a very limited 
amount of plasticity and strain at the crack tip. Therefore, the start 
of the crack extension would be stress-induced. The fundamental basis of 
fracture toughness in such a case is outlined in Fig. 11 and discussed 
below. After that, the new approach to the determination of which 

is consistent with the fundamental basis of fracture toughness, is 
outlined. 

The stress intensity factor is given by 

K ■ 0 /Jif at r « a (21) 

yy 

where a is the local normal stress at the crack tip in the plane of 

yy 

fracture and r is the distance measured from the crack tip (as shown in 

Fig. 11). Vfhen the width is large and the plasticity is limited, the 

c^yO. l//r assumed in the analysis earlier, is nearly equivalent to 

l//r stress distribution assumed in Eq. (21). If the width is small, 

the strains near the crack tip and the plastic zone is large and this 

decreases the yielding constraint at a given value of K. 

At crack initiation, a^* at r - r* (see Fig. 11) and, 

correspondingly, K - K._ - o* /2irr* where a* and r* depend on the 

It yy yy 

strength and the microstructure of the material. A material has a 
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well-defined and reproducible fracture toughness if o*^ and r* have 
clearly defined values. 

Kjg defined above depends on the magnitude of o^ at r * r* and 
this, in turn, depends on the yielding constraint. Thus, to obtain a 
reproducible value in a material, the should be measured in a 

situation where the yielding constraint is above a certain minimum value. 
The minimum value of M is probably about 3. 

The results presented in the earlier sections can be used to evolve 
a new approach to the determination of the as defined above. It 

is shown in Fig. 5 that the constraint to crack tip yielding Increases 
with the width of a CT specimen. If the width is sufficiently large, 
that is, W > Wp, a so-called plane strain fracture under high yielding 
constraint can be produced. It is also shown from a combined consider- 
ation of Figs. 4 and 5 that all points on line A have a constant value 
of M for a given a/W. At W < W^, « W and the 5Z secant 

deviation from linearity is produced only by the growth of the plastic 
zone and gives rise to line A as shown in Fig. 6. At W * V^, the 
crack growth starts and Kg ■ This is consistent with the defini- 
tion of as outlined in Fig. 11. At W > the deviation from 

linearity is produced by the growth of the plastic zone as well as the 
crack. The analysis, which evaluates the relative contribution of the 
growth of plastic zone and the crack to the total deviation from linear- 
ity and produces the line B, is shown to be essentially correct in Figs. 9 
and 10. Based on these observations, the new approach to the determination 
of Kjg is outlined below. 
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First, draw line A, with a slope of 0.25, in the versus 

W space as shown in Fig. 12. Measure Kq in a compact tension specimen 
with W > 4('^-/a = W_ and a/W ® 0.5, following a procedure identical 

to that in ASTM E399. Next, locate the measured Kg value in Fig. 12. 
The Kg value is indicated by the point marked "X". Then, draw a line B 
with slope « 0.05 through the point X. The intersection of lines A and B 
gives a Kg value equal to The approach as outlined above applies 

to the case where a Type-1 load-displacement test record is obtain' 

The operational definition of [K_^] ^ v j i 

IC 399 is based on 2% crack exten- 
sion. This is different from the K as defined in Fig. 11. In most 

xc 

instances, tKj^J^gg will be higher than the K^g as defined here. 

Table 3 compares the K__ as defined here and [K_„] for some of the 

IC 1C 3 9 9 

materials investigated by other workers [ 4 ^ 8 ]. 

Growth of the Plastic Zone and the Crack in the Three Cases 

During fracture toughness testing of specimens with width W > W^, 
the plastic zone can grow in three possible sequences: (a) the plastic 

zone grows only before the 2 % crack extension; (b) the plastic zone grows 
both before and during the 2 % crack extension; and (c) the plastic zone 
grows only after the 2% crack extension. These three possible cases 
which can be encountered during the testing are discussed below. In the 
testing situation where W < W^, Kg values obtained by the secant tech- 
nique will lie on line A. 

Case 1 . In this case, the plastic zone grows only before the 2% 
crack extension. During the 2% crack extension, the growth cf plastic 
zone is negligible. The Types II and III load-displacement test records 
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[ 1 ] which exhibit abrupt start of crack extension are usually observed In 

this case. The value of K _ can be obtained by the Intersection of 

1C 

line A and the line C which has a slope « 0.016 In the vs W plot. 

However » since the plasticity preceding crack extension may be less than 

2%, the K Is more approprl* tely determined In this case» following 
1C 

the procedure described In ASTM E399. The two procedures yield, 

though, only marginally different K values. In the limit, ig -ase 

1C 

where the p.^astlc zone grows by an amount equivalent to 2% crack extension 
prior to the abrupt physical crack extension, the value can be 

obtai> ^ d by the Intersection of lines A and C. 

Case 2 . In this case, the plastic zone grows before as well as during 
the 2% crack extension. In such a case, usually a Type-l load-displacement 
test record [ 1 ] Is obtained. The experimental data on line B In Figs. 9 
and 10 belong to this case and is obt alned by the Intersection of 

the lines A and B as shown In Fig. 12. 

Case 3 . The plastic zone does not grow before but does so only after 
the start of crack extension. Obviously, this situation Is physically 
Improbable and Is therefore a hypothetical case. However, the In 

this case should be ootalned by the Intersection of line B with the 
X axis. For a finite Kq value, such an Intersection occurs at a nega- 
tive critical width, value, and therefore has no physical significance. 

Comments on the Approach 

The approach is simple, straightforward, and is identical to 
ASTM E399 in terms of the pre .edure of measurement. 
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The procedure as outlined above can be used to determine the tough- 
ness in A compact tension specimen with crack lengths different from 
a/W » 0.5, Bowever, line B in that case will have a different slope. 

The analysis to determine lines A, B, and C has also been made for 
three-point bend and single-edge notched textstoa specimens. The gener.il 
trend of the results for the three different specimens is sisdlar; how- 
ever, the slopes of the lines are different in the three different 
specimens. 

The toughness measured here is based on the growth of the plastic 
zone equivalent to 2X crack extension. This is done for several rea- 
sons. First, it is consistcmt with E399 procedure. Second, if the 
plastic zones were larger than the 2Z equivalent crack extension, the 
domination of the l//r singularity implicitly assumed in the defini- 
tion of K may be in considerable error. Besides, a large plastic 
zone nay lower the yielding constraint too much to produce a stress- 
induced plane strain fracture. Finally, if the toughness were to be 
based t plastic zone sizes smaller than 2Z equivalent crack extension, 
an accurate determination of would be difficult because, in such -i 

case, Che corresponding or the secant values will be smaller chan 

Chat given in Table 1. It can be seen from Fig. 4 that, initially, 

K/i’y/w rises sharply with AV^/V; therefore, when the secant or the 

AV /V value is small, Che experimental values measured can be in con- 

P V 

slderablc error due to a small error in Che secant measurement. 

It should also be noted CnaC it is possible to measure in a 

specimen with W < through a plasticity corrected K value (as in 
the R-curve approach) provided the crack initiation is identified by an 
independent technique. 
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The approach to the detereination of as presented here has 

several advantages. First, It gives a value that is independent 

of the width and also of the thickness of the speciaen. The as 

defined in this approach is consistent with and R-curve approach 

since all of then are based on the start of crack growth. Finally, the 
approach recognizes that a crack can extend wider linear elastic condi- 
tions in thin and wide plates and therefore the deteraination of •4c 
in such plates is possible. Obviously » this increases the range of appli- 
cability of the approach to a wider coabination of Materials and 

configurations • 

Ia 4 >lication of the Results 

The agreeaient between the experiaental data obtained fron specisens 
with widely varying thicknesses and the predicted results in Figs* 9 and 
10 indicates that the growth of the plastic zone before and during the 
crack extension is not significantly influenced by the thickness of a 
coBpact tension specimen. On the other handg in the coopact tension 
specimen, the plasticity depends significantly on the width* The width 
dependence of plasticity, however, depends on the specimen configuration: 
for example, the analysis reported here shows that in single-edge notch 
tension specimens, the width depend^'nce of plasticity is significantly 
less than that in a compact tension specimen. In a center notch speci- 
men, where the stress gradient over the whole ligament along the direction 
of the crack is smaller, the width dependence of plasticity is expected 
to be even less. A precise, three-dimensional, elastic-plastic, finite- 
element analysis is required to understand the relative contribution of 


25 



the %rldth and the thickness of a configuration to the limited plasticity 
effects. Alternatively, the elastic and the plastic part of the thickness 
direction contraction as dependent on the thickness, width and configura- 
tion needs to be experinentally determined. 

The loading of a thin, as against a thick, structure or specimen 
often produces an out-of-plane bending, and this can obscure the eiqper- 
Imentally observed relative effects of width and thickness on the plas- 
ticity. Such bending can produce a significant amount of loading in 
the Node III at the crack tip. The growth of the plastic sone at the 
crack tip under Node III shear stress is <tuite different froci the situ- 
ation where it is absent. In fact, the use of guide plates has consid- 
erable effect on the experimental values of thin specimens. So 
li g as buckling and out-of-plane bending of the specimen or the struc- 
ture can be avoided, the width is expected to exert a relatively large 
influence on the plasticity and constraint in a given configuration. 

The width dependence of the constraint, as shown in Fig. 5, implies 
that the start of the growth of a crack in a thin but wide structure can 
occur under linear elastic condition. Thus, a R approach would be 
adequate to characterize fracture In such a situation. This is, of 
course, done at present through the R-curve approach. A J approach will 
be significant only with the width dimension U < W^. 

It may also be noted that, in the usual representation of fatigue 
crack growth, according to Paris convention, the start of stage-III 
crack growth is associated with the onset of considerable plasticity. 

The width dependence of constraint shows that the linear elastic behav- 
ior in a wider structure is extended to a higher K value. Therefore, 
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in a wider speciaen or structure, the stage-ll fatigue crack growth behavior 
will be extended to a higher AK versus da/dN reglae, whe>'e AK is the 
range of stress intensity factor and da/dN is the increment in crack 
length per cycle during fatigue loading. This has significant ramifica- 
tions with regard to the prediction of the behavior of the structure from 
the data obtained from the small specimens, particularly under spectna 
loading. 

In testing of materials, a specimen of lomr thickness, but larger 
width, may be tested to achieve a high constraint and to obtain a valid 

value. Since increases directly with B, but increases only 

with the square root of U, a thin but wide specimen will require a 
machine with smaller capacity and yet meet the constraint requirement 
for valid measurement. However, adequate care must be taken to 

meet the problems of alignment and buckling encountered during tlw 
testing of wide specimens. 

SUItlARY AND C(HICLUSIONS 

1. A method of determination of is proposed and is verified 

by comparing experimental and calculated results. 

2. The method gives a value that is independent of width. 

3. The value' determined is defined on the basis of the 

start of crack extension and is therefore consistent with the R-curve 
and approaches . 

4. The approach enables one to determine R^^ in thin but wide 
plates, and therefore increases the range of applicability of the 

R - R^^ approach to a wider combination of materials and configurations. 
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TABLE 1 — Deviation from linearity (AV^/V) due to the growth 
of plastic zone equivalent to 2Z crack extension 
(Aa/a ■ 0.02) in CT speclaeas. 



AV /V for Aa /a - 0.02 
p e p 


a/W 

from Eq. (7) 

from Eq. (7) and Fig. 4 

0.4 

0.038 

0.52 

.5 

.050 

.5 

.6 

.069 

.5 


1 

Kq is defined as the stress Intensity factor corres- 
ponding to IX crack extension as measured by the secant tech- 
nique. In the absence of physical crack extension, the 2Z 
crack extension corresponds to Aa/a ” Aa /a ” 0.02. 
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In a ICT speclaen prepared froa ASTM 516-grade 70 steel ; 
» 47 Ksi. 


e/W 

*Q’ 

K^(Bax), 


K^(max)/a^/v 


Ksl/in. 

Ksi/In. 



0.402 

36 

16 

0.54 

0.24 

0.501 

35 

17.2 

0.52 

0.25 

0.507 

34 

16 

0.51 

0.24 

0.661 

37 

20 

0.55 

0.30 

0.652 

44 

27 

0.66 

0.40 

0.667 

47 

30 

0.70 

0.45 

0.570 

49 

32 

0.73 

0.48 


1 

K^(iiax) Is the maxlawm value of the stress Intensity used 
during the final stages of fatigue cracking. 
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TABLE 3-CoaparisoD of with determined by 

the Intersection of lines A and B. 


Reference 

Material 


•4c 

13] 

T1-6A1-4V- 
Plate A 

82-92 MHb”3/2 

77 

[3] 

T1-6A1-4V- 
Plate B 

— — 

69 MNm“^/2 

[S] 

2219-T851 

Al-Alloy 

35.5 Ksi*^. 

29 Ksi/in. 
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NOTE: IN THE ANALYSIS. THE ORIGIN OF THE ORDINATE IS TAKEN 
AT THE POSITION OF THE STRAIN REVERSAL POINT. ACCORDINGLY, 
r, THE DISTANCE FROM THE CRACK TIP - X2 - X. 


FIG. 1 - Strain distribution In the ligament, (W-a), of a compact 
tension specimen . 
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NOTE: m IS A MEASURE OF THE YIELDING CONSTRAINT, AND IS 
DEFINED IN EQUATION (2). AS SHOWN IN FIGURE 1,THE STRESS 

Cyy a My/r. 


FIG. 2 — The effective and the nominal plastic zone sizes In plane 
strain condition. 
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FIG. 3 — Load-displacement curve shoving the deviation from linearity 
at any given value of P and V . 


35 


1.0 I- 



FIG. 4 — The parameter K/o„v^ as a function of the ratio AV /V , 
Y p e 

a measure of the deviation from linearity due to the formation of plastic 

zone. 
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EXPERIMENTAL DATA 
O Ni.Cr-Mo STEEL (31 
• SAE 1020 STEEL (UNPUBLISHED 
VVORK BY AUTHOR] 

X ASTM516^r70 STEEL (UNPUBLISHED 
IVORK BY AUTHOR] 

O 2219-T851/AS ALLOY (6] 

□ 7475-T7351 AS ALLOY (3] 

■ 2024 AS ALLOY (3] O 


O 


COMPUTED LINE A - SLOPE = 0^ 


.08 


.12 


.16 


.20 


W. m— 


FIG. 6 — CalculaCed and experimental values of the parameter 


Kq/Oy*'^?. 
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FIG. 7 - Schenatlc of a typical R curve for W > W^. 
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LINE C. CRACK GROWTH ONLY - 
SLOPE « a064 


\ 

LINE A. PLASTIC ZONE ONLY - 
SLOPE* 1 


I 


W/W, = 100 A*/W, 


FIG. 8 — Scheaatic representation of crack growth In the absence of 
plasticity In a CT specimen with aA^ * 0.5 . 
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FIG. 10 — Experimental and predicted results In the (K /K versus 

Cj J.v> 


W/W^ plots. 
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FIG. 11 - Definition of 
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FIG. 12 — An outline of the proposed approach for the determination 

of K__ in a CT specimen with a/W * 0.5. 

10 
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